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We give a brief introduction on finite temperature phase transitions in lattice 
QCD including a discussion on the identification of first order transitions. We 
present a study of the deconfinement phase transition of one-fiavour QCD, using 
the multiboson algorithm on lattice of sizes 8^, 12^ and 16'^ X 4. For heavy quarks 
our results are characteristic of a first order phase transition which gets weaker 
as the quark mass decreases and ends at a critical value of k ~ 0.1 or in physical 
units at about 1.6 GeV. 



1. Introduction 

Twenty years ago Polyakov§ and Susskindi pointed out that, when the tem- 
perature is increased, a deconfinement phase transition from normal hadronic 
matter to quark-gluon plasma can occur. Understanding the properties of 
QCD under extreme conditions of high temperature and/or pressure has been 
ever since a challenging problem (for a recent review see for example Ref.u). 

The nature of the deconfinement phase transition has far-reaching phe- 
nomenological consequences: In Astrophysics it is important because such a 
phase transition is believed to have occurred, in the opposite direction, 10~^s 
after the Big Bang. If of first order, effects due to supercooling should be vis- 
ible today, possibly as deviations in the light element abundance from values 
obtained in the standard scenario. Determination of the order is also of im- 
portance to model builders of neutron stars and supernovae. In experiments 
with ultra-relativistic heavy-ions the aim is to create and detect the quark - 
gluon plasma. A first order phase transition is generally considered easier to 
detect and many of the proposed signatures assume the existence of a mixed 
phased. In the planned heavy ion experiments that will start next year at 
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RHIC, Brookhaven and later on at LHC, CERN, the temperature reached will 
be of the order of 600 MeV and at this temperature one is still dealing with 
a strongly interactive system. Thus lattice QCD provides the most suitable 
non-perturbative approach to study such phenomena starting directly from the 
QCD Lagrangian. In order to use lattice QCD to study these phenomena we 
need an efficient method for simulating an odd number of flavours, namely the 
two light u and d quarks and the heavier strange quark. The locaLbosonic 
algorithm originally proposed by Liischer for two degenerate flavoursEl can be 
generalized to any number of flavours iR Before treating theJ^2 + l)-flavour 
case, we have applied it first to the study of one-flavour QCD&B. In part be- 
cause of algorithmic difficulties, one-flavour QCD has been largely ignored in 
spite of its interesting properties. The continuum Nf — 1 theory contains no 
pions and is expected to have no chiral phase transitioi0, whereas the heavy 
quark regime of the theory is expected to be qualitatively similar to the three 
dimensional three-states Potts modelEJ in an external magnetic field. 

In this talk, after giving a brief overview of how one studies first-order 
phase transitions on a finite lattice taking as an example the well-studied de- 
confinement phase transition in quenched QCD, we will present results for 
heavy to moderately heavy Wilson fermions for one flavour QCD, on lattices 
of size 8^, 12^ and 16^ x 4. 

2. Overview 

The deconfinement phase transition was first studied in the case of pure gauge. 
The basic degrees of freedom here are the SU{3) link variables Ufi{n) located 
on the link leaving site n in the direction fi oi a. space-time lattice with 
spatial size and temporal size N-t, with Ng ^ Nt if possible. The temperature 
is given by T = {Nta)~^ with a the lattice spacing. The partition function in 
Euclidean space is given by 

Z = /" [dU] e-^« and Sg = ^^[1-^ Re Tr Up] (2.1) 

where Up is the product U^{n)U^{n A- fi)Uj^{n -\- i))U},{n) of 4 links around an 
elementary plaquette P, and we sum over all plaquettes P. We take periodic 
boundary conditions for the links and, for a given Nt, the temperature is 
changed by varying /? = G/g^ with g the coupling constant. Svetitsky and 
Yaffe tHi identified the global symmetry that is spontaneously broken at the 
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deconfincment phase transition as the center symmetry Z{3) of SU{3). The 
pure gauge action is invariant under Z{3). To see this let us transform all 
time-like links with a given time coordinate no, 

C/o(n,no) zL/o(n,no) (2.2) 

with z an element of Z{3). A temporal plaquette at time uq necessarily has two 
transformed links so that Up ^ Ujx{n)zUQ{n + pL)U^^{n + Q)ul{n)z'^ = Up since 
z commutes with all t/'s. While local observables are invariant, the Polyakov 
line defined as the product of temporal links 



1 ^' 

i(n) = -Tr n C^o(n,7io) , (2.3) 

"0 = 1 

transforms as L —t zL under Z(2)). Its expectation value {L) = if the symme- 
try is unbroken whereas {L) 7^ if the symmetry is spontaneously broken. If 
(L) = the free energy of a static quark, Fq — — log((i))//3, becomes infinite 
and at large distances the correlator G(r) = (L(n)L(n -|- r)) ~ exp(— /3Fgg) 
should also vanish with an exponential dependence ^ exp(— tTa^|r| A'j) with a 
the string tension and Fqq the free energy of a static quark-antiquark pairEi. 
This behaviour describes the confining phase. If the symmetry is spontaneously 
broken then {L) ^ and Fq is finite signaling quark liberation. Therefore the 
Polyakov loop is an order parameter for the Z{?i) symmetry. Constructing an 
effective action for L and using universality arguments the authors of Ref. 
predicted that the phase transition in pure^auge is first order. This has been 
confirmed by numerical lattice simulations ej. 

Dynamical quarks are included using the standard Wilson discretization 
procedure. Writing the fermion matrix D — \ ~ kM , where M is the hopping 
matrix, the Wilson action is 

1 1 ^ 

4"^^ = J2 ^HDin, n')V(ri'); 2mga = - - — ; M = ^ M^; 

nn' fi—0 

M^{n, n') = (1 - 7M)t^M(^)'5n',n+/. + (1 + 7p)C^/,("')'5n',n~A • (2-4) 

The hopping parameter n is thus related to the quark mass, ruq, and Kc is 
the value of k for the massless limit. For the fermions we use antiperiodic 
boundary conditions in the temporal direction and periodic in the spatial di- 



rections. Under the Z{3) transformation (2.2), ip{n, no)J7o(n, no)iJj{n, no-t-0) 
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symmetry breaking can weaken the first order phase transition, change it to 
second order or even make it disappear. These possibihties must be inves- 
tigated by detailed calculations as a function of the quark mass. The order 
parameter (L) is no longer zero even in the "confined" phase since we can have 
qq pair creation shielding colour. However it will still show a discontinuity at 
the phase transition and can still be used as a probe for deconfinement. In the 
partition function the Wilson fermions can be integrated out yielding 

Z = J d[U]exp{^S°J[U]) with Sf ^ SG~TT\og{l~ kM) . (2.5) 

S(? can be expanded in powers of k and for finite temperature the k^* term 
gives a contribution to of the formEl 

~~hY^ Tr Re L(n), h>0 (2.6) 

n 

The effective magnetic field, h cx k^* to leading order, induces a magneti- 
zation along the positive real axis. Since the transition at ft, = is first order 
the magnetic term is expected to weaken it and above a critical value hep the 
transition may disappear. 



Such a behaviour was observed in the three 
dimensional three-state Potts model in an 
external magnetic fieldEl h: 

S^-Y^ /3Re<zz-/i ^ Re2„ (2.7) 

<nl> n 

with Zn an element of Z{3). Whereas Fig. |l| 
may indicate the expected qualitative be- 
haviour for QCD with a dynamical quark, 
the question of existence of an end-point 
and its location can only be answered after 
a detailed calculation. 



1 St order 




Figure 1: Phase diagram found in 
the three dimensional three-states 
Potts model in a magnetic field h 



3. Local bosonic algorithm for one flavour 

The generalization of the local bosonic algorithm!! to any number of flavours 
is made easier by finding a polynomial approximation to the fermionic matrix 
itself Hfl rather than to (75!))^. This can be done by constructing a polynomial 
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of even degree n defined in the complex plane with complex conjugate roots 
Zfe such that lim„_+oo Pn{z) — ^ for any z in the domain containing the spec- 
trum of D (not including the origin) . Since the spectral radius of the hopping 
matrix M is bounded by 8 in the free case, less in the interacting one, we are 
guaranteed that the spectrum of D = 1 — kM will remain in the complex right 
half- plane for the heavy to moderately heavy quarks we simulate (k < 1/8). 
Using the property D = 75-0^75 we can write 

n/2 

detPn{D) ^Cn\{ det{D - Zk)Uet{D - Zk) (3.8) 
fe=i 

with c„ an easily computed constant^ obtaining a local representation for detD 
detD^ lim det-^PniD)= / ff d^d^^e" S^=i (3.9) 

k=l 

The algorithm is made exact with a global Metropolis test SB at the end of a 
"trajectory" . The number n of bosonic fields (j)k is chosen so that the correction 
term leads to an acceptance rate of about 2/3. 

For the local updating of the gauge and boson fields we used standard heat- 
bath and over-relaxation algorithms as in Ref.0. A trajectory is a symmetric 
combination of {2m+l) over-relaxation steps applied alternatively to the gauge 
and boson fields, preceded and followed by a heatbath on the bosons. Ergod- 
icity for the gauge fields is maintained due to their coupling to the bosonic 
fields. The roots Zk are distributed on the circle of radius 1 centered at (1,0). 
We use even-odd preconditioning to lower the number of bosonic fields needed 
for a given accuracy, and a quasi-heatbath to initialize the boson fields using 
thermalized gauge configurations from other k and f3 valuesu. 



4. Order parameters 

In section (2) we argued that the dcconfinement phase transition is expected 
to be first order for heavy quarks. Since the discontinuities which characterize 
a first order phase transition in the continuum are smoothed out on a finite 
lattice, we must rely on finite size scaling to identify a first order transition. 
Finite size scaling was first applied to identify transitions in spin systemscJ and 
later ±0 identify the order of the dcconfinement phase transition in quenched 
QCDIIJ. The observables that we consider are based on the Polyakov loop (2.3). 
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Our strategy is to vary /3 at a given quark mass for our three spatial lattice 
sizes, and look for the following signals: 

Coexistence of the two phases: A distinctive feature of a first order transition 
is phase coexistence. On a finite lattice we can look for tunneling between 
confined and deconfined phase, which occurs over a small temperature range 
around the critical temperature. We observe enough tunneling events to study 
the double peak distribution of the norm of the Polyakov loop defined as 

17=i^L(n) (4.10) 

n 

with V the spatial volume. 

Deconfinement ratio: The deconfinement ratio is defined by 

3 1 , , 

P=-,P--, (4.11) 

with p the probability for the complex Polyakov loop to be within 20 degrees 
of a Z3 axis. Therefore if the Z(3) symmetry is unbroken we find p = 1/3 and 
p = 0, while if it is broken in such a way that the Polyakov loop is distributed 
around one axis (here the real axis) , then p — 1 and p — I. The value p = 
is obtained only in the pure gauge case where the Z{3) symmetry is exact. 
Since heavy quarks break Z(S) we obtain a value of p different from zero in the 
confined phase but still look for a discontinuity across the phase transition. 
The peak value of the susceptibility: The susceptibility gives a measure of the 
fiuctuations of the Polyakov loop. We consider the behaviour of 

xL^vam^mr) (4.12) 

which diverges at criticality for a first order phase transition in the continuum. 
On a finite lattice the discontinuity is smoothed: the distribution of xl has a 
width proportional to 1/V and a peak value XzT"'' ^ This scaling behaviour 
changes for a second order transition to x^l"^^ ^ with a < 1. For a crossover 
behaviour, where no discontinuity occurs even in the thermodynamic limit 
y — > 00, X^^""^ remains constant as V increases. 

5. Results 

The parameters of our simulations are included in Table 1. From the values 
listed in the table, we find that n cx logivolume) for a fixed quark mass and 
approximately n oc l/m^ for a fixed volume, as expectedQ. 
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Our results for the observables 
which probe the order of the transi- 
tion are shown in Figs. (|^^- Tun- 
nehng between the confined and the 
deconfined phase is clearly observed 
for K = 0.1 as shown in Fig. ^ . Simi- 
lar results are obtained for k = 0.05 
whereas for k = 0.14 tunneling is 
no longer observed. In Fig. |^ we 
show the deconfinement ratio p ob- 
tained using reweighting lis for k = 
0.05, K = 0.10 and k = 0.12. We in- 
clude for comparison the pure gauge 
results as well. Across the transition 
region the slope of p increases with 
the volume for k — 0.05 signaling a 
first order transition in contrast 
to the behaviour of p for k = 0.12 indicating that k — 0.12 is already in 
the crossover region. For n — 0.10 it is not very clear and one would need a 
bigger lattice comparable to the correlation length to make a definite decision. 
The volume dependence of the peak value of xl is displayed in Fig. |[ The 
lines shown are best fits to the form V". For k = 0.05 the best fit yields 
a — 0.96(4) whereas for k = 0.14 a = and thus we have further evidence 
that for K = 0.05 the transition is first order whereas for k = 0.14 we have a 
crossover. For k = 0.10 and k — 0.12 the situation is less clear. The small 
value of a — 0.22(3) at k ~ 0.12 as well as the absence of tunneling leads us 
to conclude that k — 0.12 is in the crossover region. For k = 0.10 tunneling is 
still observed but a — 0.56(3) so that we may conclude that we are near the 
end point of the first order phase transition. 

The critical values /3(k), given in Table 1 are obtained from the position 
of the peak value of xl- For k = 0.05 and k = 0.1 the value of f3 where the 
deconfinement ratios for various volumes cross yields a value consistent with 
the one obtained from the position of Xl^° • Not surprisingly, the shift /3c('*) — 
(3c{k = 0) from the pure gauge critical coupling is approximately half that of 
the two-flavour case. Taking the end-point value of k, K^p ~ 0.1 with 20% 
uncertainty, we can approximately map to physical units: Using the tadpole- 
improvement property Kc{P){plaq)^/'^ « 1/8 we obtain niqa ~ 1.8, with 



/c — O. 10 





sweeps 




sweeps 





Figure 2: Time history of |f2| for k = 0.10 for 
Ns = 8 (upper) Ns = 12 (middle) and = 16 
(lower) 
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(4a) ^ ^ 22QMeV from the dcconfincment temperature. This gives ~ 
1.6 Gel/ at the end-point. 



K 


8^ X 4 12^ X 4 16^ X 4 






n/acc Ksw 


n/acc Ksw 


n/acc Ksw 




0.05 


8/0.78 18 


12/0.74 20 


24/0.83 20 


5.690(2) 


0.10 


16/0.67 45 


24/0.63 50 


32/0.67 37 


5.66(1) 


0.12 


24/0.74 55 


32/0.67 40 


40/0.69 12 


5.63(1) 


0.14 


32/0.77 60 


40/0.70 37 


50/0.67 12 


5.59(1) 



Table 1: We give the k values, the number of bosonic fields n and the average accep- 
tance ace for the three volumes studied. Ksw (in kilo sweeps) Js the total number of 
thermalized configurations used in the reweighting procedure E3. In the last column 
the critical /3 values obtained from the analysis of the results are listed. 



6. Summary and Conclusions 

In this work we have shown that the multiboson algorithm is well-suited for 
the study of one-flavour QCD for moderately heavy Wilson quarks. Using an 
exact algorithm for one flavour we have carried out a detailed finite size scaling 
analysis and determined the critical value of (3 for k up to 0.14. 

We have presented conclusive evidence that the deconfinement phase tran- 
sition is first order for heavy quarks, with the critical line ending at about 
K = 0.1, i.e. mq ~ 1.6 GeV. This is in line with phenomenological expec- 
tations. The pure gauge deconfinement transition is fairly weak, with a cor- 
relation length 0{a, few cr^^/^). This is the minimmn system size necessary 
to observe the deconfinement transition. Dynamical quarks introduce a new 
length scale, the distance where the string breaks, 0{2mq/a). Confinement 
can only be observed up to this distance. When the quark mass is lowered 
sufficiently that the second length-scale is similar to (or smaller than) the first, 
one cannot tell if the system is confined or deconfined, and the transition is 
replaced by a crossover. This occurs for 

mq - Oiafew^/2) i.e. mq - 0{l)GeV (6.13) 
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Figure 3: The Deconfinement ratio for 
pure gauge and for k = 0.05, k = 0.10 
and K = 0.12 for three lattice sizes 



Figure 4: The volume dependence of the 
peaJj of the susceptibility for k = 0.05, k = 
0.10, K = 0.12 and k = 0.14. 
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